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We performed the first scanning tunneling spectroscopy measurements on the pyrochlore super- 
conductor KOS2O6 (T c = 9.6 K) in both zero magnetic field and the vortex state at several temper- 
atures above 1.95 K. This material presents atomically flat surfaces, yielding spatially homogeneous 
spectra which reveal fully-gapped superconductivity with a gap anisotropy of 30%. Measurements 
performed at fields of 2 and 6 T display a hexagonal Abrikosov flux line lattice. From the shape of 
the vortex cores, we extract a coherence length of 31-40 A, in agreement with the value derived from 
the upper critical field H C 2- We observe a reduction in size of the vortex cores (and hence the coher- 
ence length) with increasing field which is consistent with the unexpectedly high and unsaturated 
upper critical field reported. 

PACS numbers: 74.70.Dd, 74.50.+r, 74.25.Qt 



The discovery of superconductivity in the /3-pyrochlore 
osmate compounds AOS2O6 (A = K, Rb, Cs) [1] has high- 
lighted the question of the origin of superconductivity 
in classes of materials which possess geometrical frustra- 
tion [1, Q . Interest has been predominantly focused on 
the highest-T c compound KOS2O6 which presents many 
striking characteristics. In particular, the absence of in- 
version symmetry in its crystal structure [4] raises the 
question of its Cooper pair symmetry and the possibility 
of spin singlet-triplet mixing 0, 0] . 

The pyrochlore osmate compound KOS2O6 displays a 
critical temperature T c = 9.6 K, the largest in its class of 
materials (CSOS2O6 and RbOs206 which differ only by 
the nature of the alkali ion have T c s of 3.3 and 6.3 K re- 
spectively). Although band structure calculations show 
that the K ion does not influence the density of states 
(DOS) at the Fermi level @, Q, it seems to affect sev- 
eral key properties [§]. In particular, the first order 
phase transition revealed by specific heat measurements 
in magnetic fields at the temperature T p w 7.5 K has 
been ascribed to a "freezing" of its rattling motion [Io| . 
The negative curvature of the resistivity as a function 
of temperature also indicates a large electron-phonon 
scattering [ll[. Specific heat measurements [12] sug- 
gest the coexistence of strong electron correlations and 
strong electron-phonon coupling, two generally antago- 
nistic phenomena with respect to the superconducting 
pairing symmetry. The nature of the symmetry remains 
a controversial subject in the literature. NMR [l3[ and 
/iSR data suggest anisotropic gap functions with 
nodes whereas thermal conductivity experiments [l5[ fa- 
vor a fully-gapped state. 

The peculiar behavior of KOS2O6 is demonstrated 
by its upper critical magnetic field H C 2, whose tem- 
perature dependence is linear down to sub-Kelvin tem- 



peratures and whose amplitude is above the Clogston 
limit [16]. One possible interpretation is the occur- 
rence of spin-triplet superconductivity driven by spin- 
orbit coupling [5|, |6[ . Alternatively, it has also been sug- 
gested that this behavior can be explained by the peculiar 
topology of the Fermi surface (FS) sheets of KOS2O6, 
assuming that superconductivity occurs mainly on the 
closed sheet [Hj. 

The understanding of the physics of this compound 
would greatly benefit from a detailed knowledge of the 
local density of states (LDOS). Scanning Tunneling Spec- 
troscopy (STS) is an ideal tool for this, particularly 
since it allows one to map the vortices in real space and 
also access the normal state below T c by probing their 
cores In this Letter we present a detailed 

STS study of KOS2O6 single crystals, including the first 
vortex imaging in this material. 

The KOS2O6 single crystals were grown from Os and 
KO2 in oxygen- filled quartz ampoules. Their dimensions 
are around 0.3 x 0.3 x 0.3 mm 3 . The details of their 
chemical properties as well as their growth conditions can 
be found in Ref. [4. AC susceptibility measurements show 
a very sharp superconducting transition (AT C = 0.35 K). 
Our measurements are carried out using a home-built low 
temperature scanning tunneling microscope featuring a 
compact nanopositioning stage [2l[ to target the small- 
sized crystals. Electrochemically etched iridium tips are 
used for STS measurements on as-grown single crystal 
surfaces and the differential conductivity was measured 
using a standard AC lock-in technique. 

The surface topography of as-grown samples (Fig. [3J0 
reveals atomically flat regions speckled with small corru- 
gated islands a few Angstroms high whose spectroscopic 
characteristics are noisy and not superconducting (thus 
restraining our field of view for spectroscopic imaging). 
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FIG. 1: (a) Large-scale topography of KOs 2 6 (T = 2 K, 
Rt — 60 MO); the box shows the measurement area for the 
vortex maps, (b) Spectroscopic trace along a 100 A path 
taken on an atomically flat region with one spectrum every 
1 A. The spectra show raw data offset vertically for clarity 
(T = 2 K, R t = 20 MQ). 



The large flat regions display highly homogeneous super- 
conducting spectra (Fig.QJ)), which were perfectly repro- 
ducible over the timescale of our experiments (4 months). 
We have checked that the spectra obtained by varying 
the tunnel resistance R t all collapse onto a single curve, 
thus confirming true vacuum tunneling conditions. We 
have also verified that the numerical derivative of the 
tunnel current with respect to the voltage gives the same 
spectroscopic signature as the dl/dV lock-in signal. We 
stress that all measurements presented in this paper are 
raw data. 

The lack of inversion symmetry in this compound to- 
gether with several experimental findings raises the ques- 
tion of the symmetry of the gap function. In order to 
clarify this point, we have fitted our data to several sym- 
metry models, focusing on the question of the presence 
or absence of nodes and the amplitude of any possible 
gap anisotropy. We therefore considered three scenarii 
with an approximate angular dependence of the gap, i.e. 
an isotropic s-wave (Ao), a d-wave (A cos 20) with nodes 
and an "anisotropic" s-wave (Ao + Asina^) which has 
the same angular dependence as the s-p-wave singlet- 
triplet mixed state [6[. We do not take the real topol- 
ogy of the FS 0] into account, since it comprises two 
3D Fermi sheets and is hence unlikely to have any sig- 
nificant effect on the gap structure. For an anisotropic 
gap, A(0), the quasiparticle DOS is given by N(w) ex 
|Re[((cj + zr)/vV + zr) 2 - |A(0)| 2 ) ]| where T is a phe- 
nomenological scattering rate. In addition, we included 
broadenings due to the experimental temperature and 
the lock-in in our fits. The results are presented in Fig.O 
The d-wave model can be rejected at this stage since its 
zero bias conductance (ZBC) is larger than in experi- 
ment (increasing T in the model can only increase the 
ZBC). The differences between symmetries appear much 
more clearly in the second derivative spectrum (d 2 I/dV 2 , 
Fig.[2]i) which is not surprising as it emphasizes the varia- 
tions of the DOS on a small energy scale and is very sensi- 
tive to the model parameters (in contrast with the dl/dV 
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FIG. 2: Experimental and theoretical tunneling spectra, 
(a) Normalized dl/dV spectra at different temperatures from 
1.95 to 10 K (spectra are offset vertically for clarity), (b) Pa- 
rameters for the different theoretical models, (c) Comparison 
of the experimental spectrum at low temperature and low en- 
ergy with the different theoretical models; the color codes are 
explained in (d). (d) Same as (c) for the second derivative 
d 2 I/dV 2 . 



curve). The best fit is clearly given by the "anisotropic" 
s-wave model with an anisotropy of around 30%. With 
respect to the singlet-triplet mixed state, we note that 
we do not see any evidence in our data for a second co- 
herence peak arising from spin-orbit splitting. Since the 
3D nature of both sheets implies that tunneling takes 
place in both of them, the absence of a second peak also 
rules out the possibility of two different isotropic gaps 
on separate FS sheets. Our results would however be 
compatible with multiband superconductivity with two 
(overlapping) anisotropic gaps. Finally, we see no signa- 
ture of a normal- normal tunneling channel in our junc- 
tion, suggesting that all electrons involved in the tunnel- 
ing process come from the superconducting condensate. 

To investigate the temperature evolution of the quasi- 
particle DOS, we acquired tunneling conductance spec- 
tra at different temperatures between 1.95 K and 10 K 
(Fig. [2^l). The closure of the gap at the bulk T c shows 
that we are probing the bulk properties of KOS2O6. This 
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FIG. 3: Spectroscopic traces at T = 2 K across vortices for 
a field of 2 T (a) and 6 T (b). The spectra at the vortex 
centers are highlighted in red. The spatial variation of the 
conductance is shown in the corresponding insets. 



is further confirmed by the fact that similar spectra were 
also obtained on freshly cleaved surfaces. The totally 
flat conductance spectra at higher temperature show no 
support for a pseudogap in the DOS above T c , imply- 
ing that the steep decrease in the 1/(T\T) curve around 
16 K in NMR data [l3[ must have a different origin. The 
spectra taken between 6 and 9 K (not shown) were very 
noisy. This could be explained by the proximity to the 
first order transition at T p ~ 7.5 K fiol ] . 

The BCS coupling ratio 2A max //cBT' c inferred from 
our measured gaps and critical temperature is about 3.6 
for the anisotropic s-wave case, a value slightly smaller 
than that reported from specific heat measurements [l2[ • 
However, we stress that STS is a direct probe of the su- 
perconducting gap. Our findings therefore lead us to the 
conclusion that KOS2O6 is fully gapped with a significant 
anisotropy of around 30%. 

We now focus on measurements performed in an ap- 
plied magnetic field. In the vortex cores whose radial size 
is roughly given by the coherence length £, superconduc- 
tivity is suppressed leading to a drastic change in the 
LDOS which can be measured by STM. Our measure- 
ments were performed for two fields, 2 and 6 T, over the 
particularly flat region of about 60 x 60 nm 2 (Fig. []Ji). 
Each measurement was taken at 2 K with a typical ac- 
quisition time of 40 hours. 

The results are presented in Figs [3] and [4j The vor- 
tex maps (insets of Fig. [3] and Figs [4^ and SJd) show 
the ZBC normalized to the conductance at 6 meV. Fig. [3] 
displays the spectra taken along traces passing through 
vortex cores for each of the two fields considered. The 
suppression of superconductivity and its effect on the 
conductance in a vortex core can clearly be seen. The 
vortex maps show a roughly hexagonal vortex lattice 
with vortex spacings d = 352 ± 17 A and 216 ± 21 A 

at 2 and 6 T respectively, in agreement with the spacings 
1/2 

d = (2$ /i?v / 3) expected for an Abrikosov hexago- 
nal lattice [H, i.e. 345 A and 199 A. We ascribe the 



variations in the core shapes and the deviation from a 
perfectly hexagonal lattice to vortex pinning. In partic- 
ular, the vortex identified by the arrow in Fig. 0] appears 
to be split. We attribute this to the vortex oscillating 
between two pinning centers during the measurement, a 
situation which has been seen in other compounds [231 ]. 
One should also note that the islands (surface defects) 
at the border of the measurement area (Fig. [T]) could 
influence the vortex core shapes and positions. 

In order to estimate the coherence length £ from our 
measurements, we now consider the spatial dependence of 
the ZBC. Due to the proximity of the vortices, we model 
the LDOS as a superposition of isolated vortex LDOS 
which can be expressed as N(u,r) = \u n (r)\ 2 S(co — 
E n ) + \vn(r)\ 2 8{u + E n ), where ip n (r) = (u n (r),v n (r)) 
is the wave function of the nth vortex core state and 
E n its energy. An approximate solution for the iso- 
lated vortex was given long ago [24] in which the ra- 
dial dependence of each ip n (r) consists of a rapidly os- 
cillating n-dependent Bessel function multiplied by a 
cosh _1 ^ 7r (r/£) envelope common to all states. We there- 
fore construct a phenomenological model for our 2D ZBC 
maps, <j((j = 0, r) oc N(co = 0, r), by retaining the slowly 
varying parts of the wave functions alone, i.e. 



<j(u = 0, r) = <jq + A ( cosh 



(i) 



where cfq = 0.13 is the residual normalized conductance 
at zero bias in the absence of field (Fig. [2t), A a scaling 
factor, £ the coherence length and the sum runs over all 
the vortices with positions T{ in the map. Using (pQ), we 
fitted Vi and £ over the entire map for each field, thus 
considering all imaged vortices to determine £. 

The results from the 2D fits are presented in Fig. [4j^ 
and d in map format and along traces selected to pass 
through vortex cores in Fig. 0^ and f. The traces help 
to visualize the spatial extent of the vortices and assess 
the (extremely high) quality of the 2D fits. We first ob- 
serve that the normalized ZBC between the vortices is 
slightly enhanced at H = 2 T but increases strongly at 
H = 6 T with respect to the value at zero-field (Fig. [St), 
indicating a significant core overlap. From our data taken 
at T = 2 K, we obtain £ = 35 ± 3 A and 45 ± 7 A at 
H = 6 and 2 T respectively (the uncertainties are esti- 
mated from the spread of the results obtained on several 
maps: two for 6 T and three for 2 T). Using Ginzburg- 
Landau theory, we extrapolate the corresponding T = 
values as £ = 31 ± 3 and 40 ± 6 A respectively, consistent 
with the value derived from H C 2. Furthermore, our re- 
sults indicate that the vortex size decreases with increas- 
ing field and, although at the limit of the error bars, we 
believe this trend to be genuine. In addition, this finding 
is consistent with the abnormally large H C 2'. if the vor- 
tices become smaller as the field increases, the material 
can accommodate more vortices before the breakdown 
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FIG. 4: (a), (b) Experimental ZBC maps (T = 2 K) normal- 
ized to the background conductance at 2 and 6 T respectively 
with corresponding fits (c), (d); large values (red) correspond 
to normal regions (i.e. vortex cores) and low values (blue) to 
superconducting (gapped) regions, (e), (f) Experimental ZBC 
profiles across vortex centers together with the corresponding 
profiles from the 2D fits (red lines). 



of superconductivity, leading to a higher upper critical 
field. This correlates with the observed temperature de- 
pendence of the upper critical field. 

We find that the spectra at the vortex centers are fiat 
for both fields (Fig. [3]), showing the presence of localized 
quasiparticle states in the vortex cores. However, our 
spectra show no excess spectral weight at or close to zero 
bias and thus no ZBCP which is the generally expected 
signature of vortex core states. The absence of a ZBCP 
is at first glance striking considering the large mean free 
path £ « 200 nm > f in KOs 2 6 pi]. In fact, this ab- 
sence is common to many non-cuprate superconductors, 
the only known exceptions being 2H-NbSe2 ElL 0, 
and YNi 2 B 2 C " 



27]. 



Although no definitive theory cur- 
rently exists to explain such an absence, a possible ex- 
planation assumes that the scattering rate is strongly en- 
hanced in the vortex cores. This interpretation is sup- 



ported by our numerical solutions of the Bogoliubov-de 
Gennes equations for a single vortex with an r-dependent 
scattering rate T. Furthermore, these simulations show a 
radial dependence of the LDOS which is fully consistent 
with (pQ). 

In conclusion, we have presented the first scanning tun- 
neling spectroscopic measurements on superconducting 
KOS2O6. The fitted spectra demonstrate that KOS2O6 
is a fully-gapped superconductor with an anisotropy of 
around 30%, possibly resulting from a s-p singlet-triplet 
mixed state allowed by the lack of inversion symme- 
try. We have imaged hexagonal vortex lattices matching 
Abrikosov's prediction for 2 and 6 T fields. Using Caroli- 
de Gennes-Matricon theory we extract a field-dependent 
coherence length of 31-40 A, in good agreement with the 
thermodynamic estimate from H C 2. The absence of a zero 
bias conductance peak, the apparent field dependence of 
£ and the precise radial dependence of the LDOS all call 
for deeper exploration. 
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